Effects of carrier concentration on the superfluid density of high-T c cuprates 
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The absolute values and temperature, T, dependence of the in-plane magnetic penetration depth, 
Xab, of La2-xSra;Cu04 and HgBa2Cu04+,5 have been measured as a function of carrier concentration. 
We find that the superfluid density, p s , changes substantially and systematically with doping. The 
values of p s (0) are closely linked to the available low energy spectral weight as determined by the 
electronic entropy just above T c and the initial slope of ps(T)/ p s (0) increases rapidly with carrier 
concentration. The results are discussed in the context of a possible relationship between p 3 and 
the normal-state (or pseudo) energy gap. 
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Superconductivity arises from the binding of electrons 
into Cooper pairs thereby forming a superfluid with a 
superconducting energy gap, A, in the single-particle ex- 
citation spectrum. In high-temperature superconductors 
(HTS) A has essentially d x 2_ y 2 symmetry in k-space with 
Afc = Aqcos(2(/)) [1], where cj> = arctan(k y / k x ) and Ao 
is the superconducting gap amplitude which will in gen- 
eral be 4> dependent. Changes in carrier concentration 
affect the superconducting [2-8] and normal state (NS) 
[4-7] properties of HTS and there is evidence [3-7] that 
in addition to A*, there is a normal state (or pseudo) gap, 
Ajv, in the NS energy excitation spectrum in under- and 
optimally doped samples which increases with decreasing 
doping. The maximum gap amplitude shows little vari- 
ation with underdoping even though T c is reduced [3-8] , 
in disagreement with the standard mean-field Bardeen- 
Cooper-Schrieffer (BCS) theory. This unusual behaviour 
is probably linked to the presence of A^v [7]. However, 
fundamental problems such as the origin of Ajy and its 
possible effect on the superfluid density, p s , have not been 
clearly resolved experimentally or theoretically. 

The physical quantity most directly associated with 
p s is the magnetic penetration depth A. Appropriate 
systems to investigate p s as a function of doping are 
La2- 2: Sr 2 ;Cu04 and HgBa2Cu04+,5. Both have a sim- 
ple crystal structure with one Cu02 plane per unit cell, 
can have their carrier concentration controlled, and there 
is experimental evidence suggesting the presence of Ajy 
which closes with doping [3-5,9]. Here we report in-plane 
penetration depth, A a j,, measurements for high-quality 
La 2 _ a; Sr s Cu04 (LSCO) with x = 0.10, 0.15, 0.20, 0.22, 
0.24 measured by the ac-susceptibility (acs) and muon 
spin relaxation (/xSR) techniques and for HgBa2Cu04+s 
(Hg-1201) with 5 = 0.10, 0.37 measured only by pSR. We 
find systematic changes in p s with carrier concentration 



and a correlation with Aw- 

Single-phase polycrystalline samples of LSCO were 
prepared in Cambridge using solid-state reaction proce- 
dures. No other phases were detected by powder x-ray 
diffraction and the phase purity is thought to be bet- 
ter than 1%. Lattice parameters were in good agreement 
with published work [10] . High field magnetic susceptibil- 
ity measurements showed no signatures of excess param- 
agnetic centres. The measured T c 's are 30, 37.7, 36, 27.5 
and 20.3 K for x = 0.10, 0.15, 0.20, 0.22 and 0.24, re- 
spectively. These values are also in very good agreement 
with previous measurements [10]. p,SR experiments as a 
function of T were performed on the same powders for 
x = 0.10 and 0.15. Although unoriented powders can be 
used to determine X a b by /xSR [2], the acs technique re- 
quires the powders to be magnetically aligned [11]. To 
eliminate grain agglomerates, powders were ball-milled 
in ethanol and dried after adding a defloculant. Scan- 
ning electron microscopy confirmed the absence of grain 
boundaries and showed that the average grain diameter 
was ~ 5/ito. The powders were mixed with a 5 min curing 
epoxy and aligned in a static field of 12T at room tem- 
perature. Debye-Scherrer x-ray scans showed that ~ 90% 
of the grains were aligned within ~ 2.0°. Low-field sus- 
ceptibility measurements were performed at an ac-field 
H ac = 1 G rms (parallel to the c-axis) and a frequency 
f = 333 Hz down to 1.2K. Details of the application of 
London's model for deriving A from the measured sus- 
ceptibility can be found in an earlier publication [11]. 
Transverse-field-cooled ^iSR experiments were performed 
at 400 Gauss in the ISIS, Rutherford- Appleton Labora- 
tory. The field produced a flux-line lattice whose field 
distribution was probed by muons. The depolarisation 
rate, <r(T), of the initial muon spin is proportional to 
Xab~ 2 {T) [2,12]. Checks were made to ensure that the 
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values of A a fc obtained were independent of the applied 
field. The Hg-1201 [5 = 0.10 (T c = 60K) and 0.37 (T c 
= 35K)] samples were prepared in Houston by the con- 
trolled solid- vapour reaction technique [13]. 

The values of A a b(0) derived from the acs data for 
LSCO are 0.28, 0.26, 0.197, 0.193, 0.194 pm for x = 
p = 0.10, 0.15, 0.20, 0.22 and 0.24, respectively. Here 
p is the hole content per planar copper atom. The esti- 
mated error for A a /,(0) obtained by the acs technique is 
±15% and within this uncertainty the A a &(0) values agree 
with the pSR measurements. We thus find that A~ 6 (0) 
is strongly suppressed on the underdoped side, including 
optimal doping, but there is no suppression with increas- 
ing overdoping (up to p = 0.24) in contrast to reports 
for Tl 2 Ba 2 Cu06+5 [14]. Values of A a h(0) as obtained by 
^SR for Hg-1201 arc 0.194 and 0.148 p,m for 5 = 0.10 
and 0.37, respectively. We note that 5 = 0.10 and 0.37 in 
Hg-1201 correspond to p = 0.075 and 0.22, respectively 
[15]. 

The T-dependence of X a b for LSCO is shown in Fig. 
1(a) as a plot of [l/A a6 (T)] 2 oc p s (T). Data for x = 0.10 
and 0.15 obtained by /xSR are also included for compari- 
son. Overall there is good agreement between the results 
from the two techniques. From the acs data we find that 
the initial linear term in X a b(T), characteristic of a clean 
d-wave superconductor, persists up to the highest doping 
measured (x = 0.24) in agreement with electronic specific 
heat studies on polycrystalline LSCO samples from the 
same batch as those studied here [16]. Figure 1(b) de- 
picts data for Hg-1201 powders measured only by pSR, 
including data from Ref. [17] for a Hg-1201 sample (also 
from Houston) with 5 = 0.154 (p = 0.17). As in LSCO, 
we observe a change in the shape of <r(T) oc [l/A a f,(T)] 2 of 
Hg-1201 with doping. Namely, in the underdoped region 
[1/A a fc(r)] 2 shows a more pronounced curvature. Taking 
the slope of the low-T linear term to be proportional to 
P s (0)/Aq the observed trend of [l/A a £,(T)] 2 with p would 
imply that Ao remains approximately constant in the un- 
derdoped region and decreases rapidly with overdoping. 

Figure 2 shows a comparison of the present results for 
LSCO with specific heat data taken on the same samples 
[16] where A^ was observed for x = p < 0.19. In the 
inset we observe a good correlation between [1/A a fc(0)] 2 
and [S/T(T C ) - S/T(2K)] both plotted versus Sr con- 
tent x, where S(T) is the electronic entropy obtained by 
integrating the electronic specific heat coefficient 7(T) 
from to T. The quantity [S/T(T C ) - S/T(2K)] is a 
measure of the energy-dependent NS electronic density of 
states (DOS), <? n (E), averaged over ±2fcgT c around the 
Fermi energy Ep. The effect of an energy -dependent 
DOS on the London penetration depth Al or p s (0) is 
not usually considered in standard theory which implic- 
itly assumes a constant DOS and a parabolic E(fc) dis- 
persion relation. It has been argued elsewhere [5] that 
p s (0) — 47T 2 (v x 2 g n (E)) I 'e 2 where the average is over an 
(anisotropic) energy shell Ep ± Aq. Note that this result 



agrees with the standard expression for the NS conduc- 
tivity and the usual relation between Al(0) and the real 
part of the frequency-dependent electronic conductivity 
in the normal and superconducting states ai n (u>) and 
a-f(w), respectively. Namely, Al(0) is determined by the 
area under the [<7i™(a>) — <ti s (ui)] curve in the frequency 
range < hw/2n < 2Ao. Thus the inset to Fig. 2 sug- 
gests that the strong decrease of p s (0) with x from x 
= 0.20 to 0.10 is related to the suppression of spectral 
weight with energy range Ep ± A due to the presence 
of Ajv- 

The main panel in Fig. 2 shows a correlation be- 
tween the doping dependence of the initial linear terms 
of A a h(T) and the low-T specific heat coefficient 7, both 
quantities being related to the number of excited quasi- 
particles n e (T). For low values of x, n e (T=10) is much 
smaller than expected from the T c value and this proba- 
bly implies that the average value of A (</>) is significantly 
larger than T c . The rapid rise above x = 0.20 may arise 
from the combined effects of the closure of An at x = 
0.19 [4,16], the decreasing T c plus the fact that for LSCO 
there is significant pile up of states near Ep in the over- 
doped region 0.20 < x < 0.35 [16]. 

In Fig. 3(a) we present the LSCO acs data as 
[A a fc(0)/A a h(T)] 2 versus T/T c and compare the data with 
the mean-field calculation for a d-wave weak-coupled 
BCS superconductor with a cylindrical Fermi surface 
(FS) which gives A /T c ~ 2.14 [18]. There appears to 
be a systematic deviation of the data from the weak- 
coupling T-dependence with a greater (weaker) curvature 
on the underdoped (overdoped) side. We note however, 
that particularly in the overdoped samples there is pos- 
itive curvature near T c which may arise from a small 
amount of doping inhomogeneity giving a distribution of 
T c values in this region where dT c /dp is maximal [10]. 
The effect of this is to rescale the curves with a slightly 
lower value of T c . We have modelled p s {T) using the 
d-wave T-dependence and a normal distribution of T c 
values with standard deviation of 3%, 5% and 9% for x 
= 0.20, 0.22 and 0.24, respectively. These corrections, 
plotted in Fig. 3(b), bring the curves for x = 0.20 and 
0.22 into good agreement with weak-coupling BCS with 
Ao/Tc ~ 2.14. Similar deductions, as to the magnitude 
and p-dependence of Ao/T c , were made in the specific 
heat studies on these overdoped samples [16]. However, 
the x = 0.24 sample still shows significant deviations that 
possibly reflect changes in the electronic structure. This 
would not be surprising given the changes in the FS with 
the rapid crossover from hole-like to electron-like states 
near x = 0.27 [19]. We note that the data for x = 0.24 
is in excellent agreement with a weak-coupling d-wave 
calculation for a rectangular FS [20] . 

In contrast to the overdoped samples the optimal and 
underdoped samples [Fig. 3(a)], both possessing very 
small rounding near T c , diverge significantly from the 
weak-coupling curve and in the opposite direction. Wc 
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note that accounting for inhomogeneities in these sam- 
ples will, if anything, move the curves even further from 
the weak-coupling BCS fit. 

A central conclusion of the present work is that there 
is a crossover in both p s (0) and p s (T) near p = 0.20. 
Such behaviour is characteristic of many other NS and 
superconducting properties which have been interpreted 
in terms of the presence of A at in the underdoped region. 
The rate of depression of T c due to impurity scattering 
(oc I/7 at T c ) remains constant across the overdoped 
region then rises sharply with the opening of A at, be- 
ginning in the lightly overdoped region at p ~ 0.19 [21]. 
Boebinger and coworkers [22] using intense pulsed mag- 
netic fields observe a crossover from insulating to metallic 
behaviour at T = occurring near p = 0.18 and angle- 
resolved photoemission studies show the development, in 
the overdoped region, of a full NS Fermi surface [23]. In 
this region the resistivity coefficient, [p{T) — p(0)]/T, ex- 
hibits a low-T suppression due to the opening of Ajv [24] . 
These considerations provide a compelling motivation for 
interpreting our penetration depth data within a similar 
scenario. 

The proper means of incorporating the pseudogap ef- 
fects within a realistic model, and indeed the very nature 
of the pseudogap is a matter of current debate. How- 
ever, a key characteristic of A^r is the loss of NS spectral 
weight near E^. The loss of spectral weight can cause, as 
discussed above, both a strong reduction in p s (0) and, in 
a simple model, enhanced curvature in p s (T) / p s (0) above 
the BCS weak-coupling d-wave T-dependence [25], the 
very features we observe for the optimal and underdoped 
samples. 

We note that our data are in reasonable agreement 
with earlier reports for slightly underdoped grain-aligned 
HgBa2Ca 2 Cu308+,5 [11,26] and single crystal LSCO with 
x = 0.15 [27]. In contrast to the strong p dependence we 
have found in [A afc (0)/A Qb (T)] 2 for LSCO and Hg-1201, 
studies in YBCO [28,29] reported that [A a6 (0)/A a6 (T)] 2 
scaled approximately with 1/T C for various dopings, 
at all temperatures. However, systematic changes in 
[A Q h(0)/A a b(T)] 2 with p were noted at the time [29] al- 
though these were too small to allow further analysis. 
This may simply be due to the fact that the YBCO sam- 
ples were not as heavily underdoped as the x = 0.10 
LSCO sample. We also note that YBCO is complicated 
by a mixed s+d order parameter [1,30] and the effect of 
the Cu-0 chains on the total p s [31,32]. 

In summary, using the acs and /iSR techniques we have 
obtained consistent and systematic results on the effects 
of carrier concentration on p s of monolayer cuprates. In 
the overdoped region we find a more or less constant value 
of p s (0) (up to p = 0.24), and p s (T)/p s (0) is in rea- 
sonably good agreement with the weak-coupling d-wave 
T-dependence. In the optimal and underdoped regions 
p s (0) is rapidly suppressed and above 0.1T C there is a 
marked departure of p s (T) / p s (0) from the weak-coupling 



curve. In a comparative study with available specific 
heat data we found evidence supporting a link in the be- 
haviour of p s and the normal state gap Ajv- 
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FIG. 1. (a) X a b 2 {T) obtained by the ac-susceptibility technique for grain-aligned La2- a; Sr :l ,Cu04 (LSCO) with x = p = 
0.10, 0.15, 0.20, 0.22 and 0.24. Data obtained by /iSR for unoriented LSCO (p = 0.10, 0.15) powders are also included (closed 
symbols), (b) a(T) oc A~ 2 (T), for HgBa 2 Cu0 4+ i unoriented powders with S = 0.10, 0.154 and 0.37 (p = 0.075, 0.17 and 0.22, 
respectively). The data for 5 — 0.154 is taken from Ref. [17]. 



FIG. 2. Low-T X ab - 2 (T) for Laa^Sr^CuCU (LSCO) versus x compared with the low-T specific heat coefficient 7 [16]. 
Inset: A ai T 2 (0) for LSCO compared with [S/T(Tc) - S/T(2K)] [16]. 



FIG. 3. (a) [\ a b(0) / \ a b(T)] 2 obtained by the ac-susceptibility technique for grain-aligned La2- a; Sr :c Cu04 with x = p = 0.10, 
0.15, 0.20, 0.22 and 0.24 compared with the weak-coupling BCS theory (solid line) for a d-wave superconductor [18]. (b) The 
[A a f,(0) / X a b(T)] 2 data for x = 0.20, 0.22 and 0.24 shown in panel (a) but corrected for a distribution of T c values with standard 
deviation 3%, 5% and 9% respectively (see text for details). The solid lines are the BCS d-wave T-dependence corrected for 
the respective distribution in T c 's. The curves for x = 0.22 and 0.24 are shifted vertically for clarity. 
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